Cold Collision Frequency Shift of an Optical Spectrum of a Trapped Gas 
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We develop an exact sum rule that relates the spectral shift of a trapped gas undergoing cold 
collisions to measurable quantities of the system. The method demonstrates the dependence of 
the cold collision frequency shift on the quantum degeneracy of the gas and facilitates extracting 
scattering lengths from the data. We apply the method to analyzing spectral data for magnetically 
trapped hydrogen atoms and determine the value of the 15* — 25* scattering length. 
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The broadening and shifting of spectral lines of a gas 
by collisions was among the earliest discoveries in the 
development of high precision spectroscopy jjj. The 
pressure shift, which originates in interatomic pertur- 
bations is particularly simple to interpret at low 
temperatures where the thermal de Broglie wavelength 
At = h/ '(27rmfcsT) 1 / 2 is much larger than the scatter- 
ing length a Q and the interactions arise only through 
s-wave scattering. In this cold collision regime, the fre- 
quency shift is much larger than the level broadening. 

The theory of the cold collision shift has been devel- 
oped to interpret hyperfine transitions in cryogenic hy- 
drogen masers and laser cooled atomic fountains . In 
this work we study the shift for optical excitation in a sys- 
tem that can be quantum degenerate, and apply the re- 
sults to data on 1S—2S two-photon excitation of trapped 
atomic hydrogen ||. 

For the case of a homogeneous sample of density n, 
and a coherent, weak excitation that couples two inner 
states of the atoms, we find 



g 2 (Ai2 - An)? 



A Q/ 3 = 4irh 2 a a p/m . (1) 



Here gi is the equal point value of the second order corre- 
lation function p], gi = g^ 2 '(r = 0), the state 1(2) is the 
ground (excited) state, and a a is the s-wave scattering 
length for a — (3 collisions. 

Equation (|l]) shows that quantum correlations in the 
system are manifest in the collision shift. For a uniform 
Bose gas in thermal equilibrium 52 = 2 — («becA 1 ) 2 @j 
where tt-bec is the density of condensed atoms. Above the 
condensation temperature, when jt-bec = 0, 32 equals 2, 
in which case Eq.(|l|) is in agreement with previous work 
At zero temperature, for a pure condensate with 
fiBEC = n, the collision shift is half of the shift for a non- 
condensed gas. Equation (Q) generalizes the result Q to 
T < Ibec and relates the spectral shift to the condensate 
fraction. 

It is quite remarkable that the factor 172 in Eq.([l|) mul- 
tiplies both A12 and An. This results from correlations 
between an excited atom and other atoms. During the 
excitation, the internal states of the atoms are rotated: 
cos9(t)\lS) + e-^Wsin0(*)|25). The angles 0(t), <f>(t) 
depend on laser power and on the atom's trajectory in 



the laser field, specific for each atom. However, for small 
excitation power, the angle 9(t) is small, and thus the 
internal states of all atoms remain nearly identical while 
the laser is on, even if the excitation field is spatially 
nonuniform. Therefore, during the excitation the atoms 
interact as identical particles. This causes the short range 
statistical correlations in the initial state to be replicated 
in the excited state of the gas, which results in the sta- 
tistical factor <?2 in the first term of Eq. (jl|) . 

The transfer of spatial correlations to the excited state 
is not limited to weak excitation. For the case of strong 
excitation, spatial correlations in the ground state will 
also be transferred to the excited state, but only pro- 
vided the excitation scheme is coherent. The difference 
between coherent and incoherent cases can be seen from 
comparing two examples, the coherent superposition of 
the ground and excited states obtained, e.g., by a ir/2 
pulse, and the incoherent mixture state resulting from 
saturating the Rabi transition. These states will both 
have equal populations in the two internal states, but 
quite different correlations. In the former case of a pure 
internal state the spatial correlation will be the same as 
for the ground state of indistinguishable particles. In the 
latter case of a mixed state the correlations will be re- 
duced. Consequently, the correlation energy of the first 
state will exceed that of the second state by the factor 

52- 

To emphasize the non-trivial character of the result 
Eq. ([j]) , let us point out that Ti Aw co u differs from the ther- 
modynamic work needed to transfer one atom from the 
state 1 to the state 2. The latter work, calculated by 
removing one atom from the sample, and then introduc- 
ing an atom in state 2 from far away, ignoring entropy, 
is given by (A12 — <?2An)n. Here Ai2« is the energy of 
interaction of the excited atom with the atoms in the 
state 1, and 02Ann is the chemical potential of a Bose 
gas. The key difference between this process and optical 
excitation, resulting in the different dependence on 02, is 
the incoherence of the state of the added atom with the 
initial state of the sample. 

To calculate the full optical spectrum shape of a 
trapped gas in the cold collision regime, other factors 
would have to be considered in addition to the effects 
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of statistical correlations. Optical coherence can be lost 
via dephasing elastic collisions, giving rise to collisional 
broadening. One would have also to take into account 
atomic motion in the trap and the effects of the inho- 
mogeneous density distribution in the sample, especially 
in the Bose-Einstein condensate. In addition, the in- 
teraction may give rise to a doublet structure of the 
spectrum 0. Altogether, these effects can lead to a 
complicated broadened spectrum with asymmetric lines 
P, fH)|]i~i] ] . However, we demonstrate below that the spec- 
trum's center of mass obeys a simple and exact sum rule 
and is insensitive to these additional effects. 

We lay out the theory of the shift by deriving a sum 
rule (Eq.(|l2"|)) that relates the center of mass of the ob- 
served spectrum to measurable experimental parameters. 
The sum rule bridges between the uniform density re- 
sult Eq. (jl|) and experimentally measured spectra. The 
sum rule accounts for all interactions between atoms oc- 
curing in the s— wave scattering channel, which includes 
the s— wave collisional broadening. It follows from the 
sum rule that collisional broadening as well as the time 
of flight broadening resulting from atomic motion in the 
trap do not contribute to the spectral shift. At the same 
time, the effects on the shift of inhomogeneity in the gas 
density and non-uniformity in the excitation field are ex- 
pressed in the sum rule Eq. ( |l2| ) in an exact and straight- 
forward way. 

We start by considering a homogeneous Bose gas and 
derive Eq.(|l|). Then for the realistic situation of a 
trapped gas sample we derive the sum rule Eq.(|l2"|), a 
generalization of Eq.(|l|). The sum rule is exact and 
general, applicable both to Doppler-free and Doppler- 
sensitive spectra. Finally we apply the sum rule to exper- 
imental data on the spectrum of cold trapped hydrogen 
to calculate the 15 — 25* scattering length for hydrogen. 
The system: To provide the context for the theory, we 
briefly describe the experimental situation. The temper- 
ature of the hydrogen is 30 — 100 fj,K, well below the cold 
collision threshold T ~ 1 K ||. The atoms are spin 
polarized and interact in the triplet channel. Calculated 
values of the 15— 15 and 15*— 25 triplet scattering lengths 
are a n = 0.0648 nm @ and a u = -2.3 nm Q. We 
neglect 25 — 25 scattering because the excitation rate 
is assumed low (in the experiment typically 10 -4 of the 
atoms are excited) so the background gas is essentially 
pure IS*. Since |ai2 3> an, collisions between 15* and 2S 
atoms dominate the shift, which is to the red. 

Each atom will be in some superposition of the ground 
state 15 and the excited state 25. In the second quanti- 
zation formalism, the atoms are described by the canon- 
ical Bose operators ipi (r) and tp2 (f) ■ The Hamiltonian is 
TC = TCq + 7ii n t , where TIq describes atoms freely moving 
in the trap, and 7ii n t is the interaction term: 

= / E i + ('')f-^ + c/w)*«W^, (2) 




Here U (r) is the trap potential (essentially the same for 
the 15 and the 25 states). 

Inelastic collisions, such as collisions in which the hy- 
pcrfine level of one or both of the colliding partners 
changes, may contribute additional shifts which are not 
accounted for in this formalism. However, these effects, 
as well as the three-body collision effects, are small in the 
experiment and can be neglected. 

The two-photon 15 — 25 spectrum consists of Doppler- 
free and Doppler-sensitive excitations. In the Doppler- 
free situation, the transition results from absorbing two 
counter-propagating photons with equal frequencies and 
zero net momentum. In the absence of interactions, the 
resonance condition is 2wi aser — ujq , where ujq corresponds 
to the resonance of a single free atom. In the Doppler- 
sensitive situation, the transition is caused by two pho- 
tons propagating in the same direction. For a free atom, 
the resonance frequency is shifted by the recoil energy: 
2fi,u;i aS er = ?iwo + (2fc) 2 /2m, where k = hco/c is photon 
momentum and m is the atom mass. 

Radiative excitation in a many particle system is de- 
scribed by adding to the Hamiltonian (H),(||) the term 

H rad = J d 3 r(^(r)e-^V 2 + M^iM + h.c.) , (4) 

where u = 2wi ascr — ujq. The two-photon excitation field 
A(r) is equal, up to a constant factor, to the square of the 
electric field. Spatial variation of A(r) in the Doppler- 
free case occurs on a scale set by the focused laser beam 
diameter, and in the Doppler-sensitive case is given by 
A{r) cos(2kr+(j)(r)), where A(r) and <fi(r) are slowly vary- 
ing functions. 

A tutorial example: Before discussing the general 
case, here we derive the mean frequency shift for the 
Doppler-free transition caused by a uniform excitation 
field A(r) = Aq, ignoring the 15 — 15 interactions 
(An =0). To that end, consider a gas of N atoms con- 
fined in a box of volume V. Since we ignore the 15 — 15 
interaction, the many body state ground state of the sys- 
tem $o is simply a symmetrized product of single particle 
states. It can be characterized by occupation numbers 
rij of the single particle plane wave states V~ 1 ^ 2 e lkjr , 
J2j n j — N. Initially, the internal state of all atoms is 
15. 

The excited state, to lowest order in the excitation, is 
given by $i = Hrad^o- We consider the norm || 4>i || 2 
and the expectation value of the interaction ($1 |Wi n t 
The ratio of these quantities gives the mean frequency 
shift. Because $0 is the product of plane wave states in 
a box, the frequency shift can be evaluated exactly. 

The norm of the excited state is given by 
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\ 2 J (&oH>t {r)ip2(r)ip2 ( r ')*f } i( r ')\^o)d 3 rd 3 r l We start with the Golden Rule formula for the absorption 



spectrum, 



(5) 



To evaluate the norm one first puts the operators ip2(r) 
and V , 2~( r ') m © m normal order by using the commu- 
tation relation [^{r), ipt( r ')] = ~ r ')- Noting that 
V^MI^o) = 0, the norm is given by 

|| $! || 2 - \A a \ 2 / ($ |^ + (r)Mr)\^o)d 3 r = \A \ 2 N . 



2n 



]T 5(hcu + E t - E f )\(f\H iad \i)\ 2 Pl , (10) 



where i), |/} are eigenstates of the Hamiltonian H = 
TCo + Hint with the energies Ei, Ef, and pi is the statis- 
tical occupation of the states 

The sum rule for the spectrum I{u>) is found by eval- 
uating the first moment: 



(G) 

To obtain the frequency shift Aw co n, we consider the ex- 
pectation value (^llWintl^i), keeping in 7i int only the 
15 — 2S interaction A12. After arranging in normal or- 
der, as in the calculation of the norm || $1 || 2 , one has 



ujX(lu 



du 
2^ 



(*i|Wint|*i> = MM' J (®om(rm(r)Mr)i>i(r)\$o 

(7) 

Evaluating the expectation value for $0 chosen as a prod- 
uct of plane wave states, one expresses Eq.(Q) in terms of 
the occupation numbers of the ground and excited states 
as 



)d 3 r 



p ( E f - mim^fPi 

^3 (<l«rad|/)(/|[W,«rad]|*)Pi 
1 

n 



Et,E } 

= ^3 ^2(i\Hi a d[H,H Ia d]\i)pi 

Ei 



(11) 



(8) 



In obtaining this result we first integrated the delta func- 
tion, then wrote the result as a matrix element of the 
commutator [H, 7i r ad] and, finally, used the completeness 
relation. 

Now we consider contributions of the different terms of 
the Hamiltonian to the sum rule. The potential energy 
operator J (tpt^i +^2 V'2) U(r)d 3 r commutes with H ra ,d, 
and thus does not contribute. There are two contribu- 
tions, first from the interaction Hamiltonian, second from 
the kinetic energy operator, denoted by Fi nt and F^ n re- 
spectively. The sum rule becomes 



The mean frequency shift is then given by the ratio of 
(0) and the norm 



HAuj, 



coll 




/duj 
uj1(uj)— = F^t + -Fkin 



(12) 



^ n t (rii + 1) 



(9) 



The formal reason for the factor 2 to appear in Eqs. (||, ||) 
and, eventually for g 2 to appear in Eq. ([!]), is the follow- 
ing. In taking the average in Eq.(^) by Wick's theorem 
jTj|, there are two essentially different ways to pair the 
operators, analogous to the Hartree and Fock contribu- 
tions to the energy. For short range interaction between 
bosons, the Hartree and Fock contributions are equal and 
as a result the frequency shift is twice as large as the 
"mean density" result. 

In the thermodynamic limit, V, N — ► 00, n = N/V con- 
stant, the second term in Eq.(^|) contributes only when 
there are states filled by a macroscopic number of par- 
ticles. For example, in thermodynamic equilibrium at 
T < T B eci the shift Eq.(|J) is Ai 2 (2n - n 2 /n), whereas in 
a non degenerate gas, at T > TbeC; the shift is 2Ai2?i. 
The sum rule: We turn now to deriving a sum rule that 
generalizes the result Eq.([l]) to non- homogeneous sam- 
ples and spatially varying excitation field (and An 7^ 0). 



For Doppler-free excitation F^ is small compared to 
F^t, whereas for Doppler-sensitive excitation it con- 
tributes the larger shift. 

First, consider the interaction 7ii n t, and calculate -Fi nt . 
After evaluating the commutator in Eq.([ll]), one follows 
the same procedure as in the above calculation of the 
norm II $1 ||. The result is 



(A12 - An) \A(r)\ 2 ^+(r)^+(r)Mr)Mr 



d 3 r 

'IF 

(13) 



where (...) means YIe- (*l---l*).Pi- The expectation value 
(: (ipf (r)ipi(r)) 2 :) = G^r), the two-particle density. 
(Here : ... : indicates canonical normal ordering.) Fi- 
nally, using the statistical factor g 2 — G2/™ 2 , the result 
is 



Fu 



(A 12 -An)|A(r)|W(r) 



d 3 r 

IF 



(14) 



Next, we calculate -Fkin, the contribution to the 
sum rule coming from the kinetic energy operator 
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-(h 2 /2m) f(ip^V 2 ipi + ■0^V 2 V'2)rf 3 r. After evaluating 
the commutator with 7i ra d, one has 

^ kin = ~£(/ ^t(r)A*(r)[V 2 ,A(r)} ^(r)^\ . 

(15) 

Integrating by parts, and writing the excitation field as 
A(r) = \ A(r)\e if> , yields 

Fkin = / (J^\VA\*n-*\A\ 2 3-Ve) ^ , (16) 

where n and j are the particle number and flux densities: 

ih 

n(r) = (i>t(r)il>i(r)), j(r) = -— W (r)Vi>i(r)) +h.c. 

(17) 

The first term in Eq.([l6|) generalizes the ordinary mo- 
mentum recoil energy shift to the trapped gas prob- 
lem EHj. The second term represents the Doppler shift 
due to possible macroscopic gas flow in the sample. 
To clarify this, consider A{r) — Aoe tpr / fi , which would 
describe Doppler-sensitive excitation. Then F^ n — 
\Ao\ 2 J(p 2 /2m — p ■ v)n d 3 r/h 3 , where v — j/n is the 
local velocity. The sensitivity of the frequency shift to 
motion within the sample, manifest in the second term 
in Eq.(|l6|), makes it possible, in principle, to detect vor- 
tices in the condensed state. 

To employ the sum rule, one needs to relate the inte- 
grated spectral power to A(r) and n(r). Repeating the 
steps that led to Eq.(^2|), one obtains 

l tot = /*M^= J \A(r)\ 2 n(r)^ (18) 

Combining Eq.(p"8]) with the sum rule Eq. (|l~2]) , one ob- 
tains an exact expression for the spectrum's "center of 
mass" to — J LoT(u;)duj/ Jj(uj)duj. 

For example, consider a uniform density sample, and 
ignore the spatial variation of the laser field A(r). 
Eq. @ gives Fin* = (A12 - An) gin 2 J \A(r)\ 2 d 3 r/h 3 . 
In our experiment -Fkin can be neglected. Simplifying 
Eq.(|i~8|) and combining it with Eq.(|l4]) yields the fre- 
quency shift Eq.(|l|). 

There are two comments concerning the generality of 
the sum rule. First, note that in deriving the sum rule 
Eq. (|i~2|) , we do not assume thermodynamic equilibrium. 
The result is exact and applies to non-equilibrum sys- 
tems for which the factor gi may differ from its equilib- 
rium value. Second, the above derivation of the sum 
rule assumes coherence of the excitation described by 
(kl|). One can see, however, that the results 
hold as well for an incoherent excitation field of the form 
A(r)e iwt+i< ^*) with a fluctuating phase <f>(t). Also, it is 



straightforward to generalize the results for the excita- 
tion field with different spatial dependence of different 
frequency components. 

Analysis of the data: To investigate the utility of the 
sum rule, we applied it to extract a value for ais_2S from 
data on the IS — 2S transition in hydrogen for a normal 
gas. An account of the experimental situation and ex- 
amples of the spectral data will be published elsewhere 
p6[ . Evaluating the integrals in the sum rule requires 
knowing the excitation field A(r), the value of g2, and 
the density n(r). The excitation field is generated by a 
Gaussian beam which is fully characterized by a single 
parameter, the beam radius, which can be accurately de- 
termined. For temperatures above TbeCj 9i = 2. At 
lower temperatures, 172 depends on the temperature. 

The density n(r) was found by measuring the peak den- 
sity and knowing the properties of the trap. The peak 
density uq was determined by exploiting the property 
that the BEC critical density is accurately described by 
the ideal gas expression: n c = 2.612(27r7t 2 /ksTm) 3 ^ 2 . 
If the system is at the transition point then measuring 
the temperature determines the density. The system was 
cooled into the condensate regime and the spectrum was 
observed as the condensate decayed. The spectrum was 
measured after the condensate had decayed for ten sec- 
onds, when the presence of a small though visible con- 
densate assured that the peak density n had its critical 
value. The contribution of the condensate to the spec- 
trum was unimportant. Thus, the only quantity required 
to apply the sum rule was the temperature. This was 
found by measuring the width of the Doppler-sensitive 
spectrum |i~7j ]. 

From the experimentally measured spectrum we found 
v = —29 ± (2) KHz. We numerically calculated the inte- 
grals on the right hand side of Eqs.(|l^),(|l~8|) and found 
that 2(Ai2 - Xu)/h = 4.4 ± (1.7) x lCT 10 Hzcm 3 , where 
the major sources of uncertainty are the temperature and 
the trap and laser geometry. From this wc determined 
the 1S — 2S scattering length to be CJ12 = — 1.6±(0.7) nm. 
This result is in reasonable agreement with the calculated 
value |l8[ , ai2 = -2.3 nm. 

An alternative approach to extracting the scattering 
length was used in Ref. p], where the value an = 
— 1.4 ± (0.3) nm was reported. The interpretation em- 
ployed a semiclassical description of the atomic motion 
and a local density description of the phase shift. The 
present method is more direct and, we believe, more reli- 
able. It can be viewed as a check on the earlier analysis, 
and a confirmation of the calculation of the dipolar decay 
constant J20| on which it depends. 

In reference |lQ ], an internally consistent description 
of the density of the condensate required assuming that 
the density shift parameter in the condensate was the 
same as in the normal gas, rather than half as large as 
expected from this analysis. This anomaly remains to be 
explained. 
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As a speculative explanation, one could consider a 
state close to the transition temperature in dynamical 
but not in thermal equilibrium, in which several low 
energy states are populated with macroscopic occupa- 
tion numbers. For example, for N particles distributed 
equally among m states, one has rij = N/m, i = 1, ...,m 
in Eq.(§). Then from Eq.(§) ?iAw coll = (2 - l/m)A 12 n, 
i.e., the shift is described by an effective #2 = 2 — 1/m. 
For a large number to of constituent states, the effective 
<72 can be arbitrarily close to 2. 

In summary, we have shown that quantum statistical 
correlations of a cold gas sample are imprinted in the 
collisional shift of the center of mass of an optical ab- 
sorption spectrum. In the cold collision regime the sum 
rule Eq.(|l^) can be applied to determine the statistical 
correlation factor from optical spectrum. The sum 
rule is valid for any gas in the cold collision regime. It 
takes into account possible inhomogeneities in the sample 
and the excitation field, and it is valid above and below 
Tbec- Also, the sum rule is valid for a non-equilibrium 
system, with 172 values possibly different from those in 
equilibrium. We have demonstrated the usefulness of the 
sum rule by using it to extract the IS — 2S scattering 
length for hydrogen from experimental data. 

It should be pointed out that our results, the frequency 
shift (|l|) and the sum rule (|T^), are only valid at small 
mixing angles of the IS and 2S states. The cold colli- 
sion shift at large angle mixing is an important problem, 
particularly for atomic clocks. The generalization of the 
results ( |l2| ) and (|l|) for such systems is an interesting 
open problem. 
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